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Abstract 

We study the weak values of a quantum observable from the point of 
view of the Wigner formalism. The main actor is here the cross- Wigner 
transform of two functions, which is in disguise the cross- ambiguity 
function familiar from radar theory and time-frequency analysis. It 
allows us to express weak values using a complex probability distribu- 
tion. We suggest that our approach seems to confirm that the weak 
value of an observable is, as conjectured by several authors, due to the 
interference of two wavef unctions, one coming from the past, and the 
other from the future. 
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1 Introduction 

We study in the present Letter the notion of weak measurement introduced 
by Aharonov and Albert [2], Bergmann, and Lebowitz in [TJ [2j [3l 01 [5j [6] 
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from the point of view of the Wigner phase space formalism. This will allow 
us to discuss the claim made by these authors that the weak value can be 
seen as the interference of two wavepackets, one going forward in time and 
the other backwards in time. 

Let us briefly recall the difference between an ideal (also called strong, or 
von Neumann) measurement, and a weak measurement. Let A be a (quan- 
tum) observable, realized as an essentially self-adjoint operator; we assume 
for simplicity that A has a eigenvalues a%,a2, ■■■ with corresponding orthog- 
onal eigenfunctions ipi,^, ■•• ■ In an ideal measurement the expectation 
value of A in a pre-selected state ip is 

if the sequence of eigenvalues lies in some interval [a m i n ,a max ] then we will 
have a mm Smax' In fact, if one performs the ideal measurement 

the outcome will always be one of the eigenvalues Xj, and the probabil- 
ity of this outcome is |Aj| 2 /||V>j|| 2 where Xj is the coefficient of tpj in the 
Fourier expansion ip = ^2jXjtpj. Moreover the system will be left in the 
state ipj after the ideal measurement yielding the value aj. The situation 
is very different for weak measurements. As is explained in Ritchie et al. 
[25j (also see Berry and Shukla |9j, Steinberg |27J), in a weak measurement 
the eigenvalues are not fully resolved and the system is left in a superposi- 
tion of the unresolved states. If an appropriate post-selection is made, this 
superposition can interfere to produce a measurement result which can be 
significantly outside the range of the eigenvalues of the observable A. The 
post-selection can then be accomplished by making an ideal measurement 
of some other observable B and selecting one particular outcome. Thus, 
the post-selected state \<j>) is an eigenstate of B which can be expressed as 
a linear combination of the eigenstates of A (we note that, conversely, an 
ideal measurement can be expressed as a convex sum of weak values: see 
Hosoya and Shikano [21|). If (<f)\ip) ^ (and if (ft, ip are square integrable) 
the weak value of A is then the complex number 

'AtL = (2) 



We will show that this weak value can be expressed in terms of the 
cross- Wigner transform 



W(<j ) ^)(x, P ) = {^ r f I e -&v ( t>*(x + \y)^(x-\y)dy (3) 
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of the pair ((f), tp) whose physical interpretation is that of an interference 
term in the Wigner distribution of the sum (f) + t/j; we mention that the 
importance of these interference terms have been emphasized and studied 
by Zurek [29J in the context of the sub-Planckian structures in phase space. 
The cross- Wigner transform is a very important object being intensively 
studied in the harmonic analysis literature and in time-frequency analysis; 
see e.g. Cohen [TT], Folland jH] Grochenig [17] . Hlawatsch and Flandrin 
[20] . Notice that W(4>, ip)(x,p) reduces to the familiar Wigner distribution 
(Hillery et al. [E], Littlejohn [22] when cp = tp. 

We will not address here the ontological debates arising around the prob- 
lem of "Elements of Reality" (see Cohen and Hiley [I2j [13] ) ; these questions 
are difficult and have led to profound philosophical controversies. 

Notation. We will work with systems having N degrees of freedom. Po- 
sition (resp. momentum) variables are denoted x = (xi, ...,xjv) (resp. p = 
(pi, -;Pn))'i they are vectors in M. N . The corresponding phase space vari- 
able is z = (x,p); it is a vector in phase space M. 2N . We will endow the phase 
space with the standard symplectic form a(z, z') = px'—p'x. When integrat- 
ing we will use, where appropriate, the volume elements dx = dx\ ■ ■ ■ dxN, 
dp = dpi • • • dppj, dz = dpdx. The unitary h- Fourier transform of a function 
ip in L 2 (R N ) is defined by 




2 The Main Result 



2.1 A complex probability distribution 

The cross- Wigner transform ([3]) satisfies the "marginal properties : 




(4) 



and 




(5) 



It follows from the equality that 




(6) 
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For (4>\tp) / Owe define 

W{^)(z) 
= (^j • ( 7 ) 

Note the conjugation relation p^^z)* = p^^{z); also p\<t>^{z) = p^{z) 
for every complex A 7^ hence the function p^ ^ only depends on the states 
\ip) and \<p). In view of Eqn. ([6]) we have 



z)dz = 1 (8) 



hence can be viewed as a complex probability distribution with respect 
to which the average of the classical observable A is calculated; also, Eqn. 
([8]) implies that 

Re p^(z)dz = 1 , / Imp^^dz = (9) 



so that Re can be viewed as a quasi-distribution, in the same way as the 
usual Wigner transform. When ip = (p then Imp^^ = and Re p^^ = Wtp. 
Observe that it immediately follows from Eqns. ([7]) and (HJ), © that the 
marginals distributions of p^ are given by 

note that anyone of these equalities allows by integrating in the conjugate 
variable to recover the normalization condition 0. 

We point out that the consideration of complex probability densities 
has per se nothing unusual; such complex probabilities have been used in 
the context of stochastic processes (see Zak [28]), signal theory (multipath 
fading channels, see Chayawan |10j ) and they also appear in the study of 
non-Hermitian quantum mechanics (see Barkay and Moiseyev [8]). 

We claim that: 

Theorem 1 Let A be a classical observable and A its Weyl quantization; 
we have 

(%tL= f A(z) PM (z)dz. (11) 



The reader familiar with the Weyl-Wigner-Moyal formalism (de Gosson 
[151 116] . Littlejohn |22j ) will have noticed that when <j> = ip formula (jlip 
reduces to the well-known relation 

(il>\A\il>) = [ A{z)W^{z)dz 
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yielding the usual expectation value (A)^ = (ip\A\ip) / . We will study 
the relative importance of these values when (f> and if) are coherent states in 
Subsection I2.3I 



2.2 Proof of Theorem 1 

To prove formula (fTT]) it is sufficient, in view of definition (J7J) of p ( f ) ^(z), to 
show that 



A\ij})= W((j),ip){z)A(z)dz. (12) 



To prove the latter we could perform a direct calculation staring from the 
right-hand side, inserting the expression ([3j) of W(4>,ip)(z) and making var- 
ious changes of variables. We prefer to give a more elegant proof which has 
some conceptual advantages. The first step consists in observing that the 
cross- Wigner transform can be expressed in terms of the Grossmann-Royer 
[151 [26] operator 

T GR (z )4>(x) = e^~ X0 U(2x - x) (13) 

(also see de Gosson [16], Chapter 9). A simple calculation shows that we 
have 

WfaMz) = {±f (f GR (z)M) (14) 
and that the Weyl quantization A of the observable is given by 

AiP(x) = (±) N [ A(z )f GR (zMx)dz . (15) 

Using the latter we have 

= [±) N I A(zoMT GR (zoWdz ; (16) 



we next observe that T GR (zo) is both unitary and involutive (i.e. Tgr(-Zo) = 
Tgr(zo) -1 ) and hence 

<0|T GR (z o )V) = {T G k(zo)- 1 ^) = (Tgr(^o)</#) (17) 
so that (I16p can be rewritten 

(4>\A\i>) = (±) N [ A(z )(f GR (z )^)dz 

= [ A(z )W(cf>,i>)(z)dz 
which was to be proven. 
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2.3 The case of coherent states 

Suppose that both wavefunctions are normalized coherent states concen- 
trated near zq = (xo,po) and — zq at time ii n , that is we choose 9 and 
tp = ip ZQ where 

= (inf 14 T(z )e-r>\ 2 , i>{ x ) = (±) N / 4 T(-z )e-^ 2 ; (18) 

where T(zq) = e~h( XoP ~ PoX ' is the Heisenberg-Weyl operator. These states 
are minimum uncertainty states (they saturate the Heisenberg inequalities 
AxjApj > | ft). A standard calculation of Gaussian integrals shows that the 
scalar product of these states is 

(8\ip) = e -rM 2 . (19) 
Let us calculate W(4>, ip). Using the translation formula (see de Gosson [16]) 
W(T(a)<j), T(P)ip)(z) = e-i x ^ (z) W((j), ip){z - \{a + &)) (20) 
where Xa/3 is the phase function defined by 

Xap(z) = l<r(z,a- P) + <r(a,P) (21) 
(a the standard symplectic form). We thus have 

W{^){z) = e L n°^W{^ti ){z) 

where a(z,zo) = pxo — pox and £q(x) = (Tth)~ N l 4: e~\ x \ 2 l h is the standard 
fiducial coherent state (Littlejohn [22]). Now, W^(£o)£o) = W&i the Wigner 
distribution of £o> which is given by 

Wto(z) = (±f e-rM 2 , \ z \* = \ x \* + \p\* (22) 
(de Gosson [EKE], Littlejohn [22]). We thus conclude that 

WtfMz) = (J H ) Af e l-(^o) e -i|,| 2 . (23) 

Using the scalar product formula (|19p we see that the complex probability 
distribution is given by 

P^) = (^) N e^ Z0) e^e-^. (24) 
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This formula shows that p a ^{z) has an oscillatory behavior which is sharply 
peaked near the origin. We notice that since 



The integral in the third line is easy to evaluate; its value is (ttK) hence 
we have the estimate 



This inequality shows that even if the observable A is small, the weak value 
can a priory take very large values provided that the phase space distance 
between both wavepackets eft, tp is large; this is in strong contrast with what 
happens for the individual states \4>) and \tp), for which lead to the estimates 



the relative phase space localization of these states does not play any role 
in these inequalities. We will shortly discuss non-trivial extensions of the 
superposition considered above in the discussion below. 

3 Discussion 

Let us apply the phase space formalism to a discussion of the situation ini- 
tially considered in [31 H] where at a time tj n an observable A is measured and 
a non-degenerate eigenvalue was found: \ip(ti n )) = \A = a) (the pre-selected 
state); similarly at a later time t^ n a measurement of another observable B 
yields \4>{t^ n )) = \B = b) (the post-selected state). Let t be some intermedi- 
ate time: t; n < t < ig n . Following the time-symmetric approach to quantum 
mechanics (see the review in [6]), at this intermediate time the system is 
described by the two wavefunctions 



the weak value (-A)^ k satisfies 






(25) 



< sup \A(z)\ , \(A)^\ < sup \A(z)\; 



(26) 
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where UK, = e %H<yt is the Schrodinger unitary evolution operator 

(H the quantum Hamiltonian) . Notice that (f>t travels backwards in time 
since t < t^ n . The situation is thus the following: at any time t' < t the 
system under consideration is in the state \ipf) = U^ t . \ijj{t\ n )) and has 
Wigner distribution Wtpt'] at any time t" > t the system is in the state 
\4>t") = Up, t \(p(tin)} and has Wigner distribution W4> t '- But at time t it is 
the superposition \ip t ) + \(p t ) of both states, and the Wigner distribution of 
this cat-like state is 

W{(j)t + i>t) = W<f> t + WiPt + 2ReW(<i> t ,iPt)- (27) 

This equality shows the abrupt emergence at time t -and only at that time!- 
of the interference term 2ReW((fit,ipt), signalling a strong interaction be- 
tween the states \tpt) and \4>t)- Such an interaction is due to the wavelike 
nature of quantum mechanics, and is absent from classical mechanics. The 
appearance of interference terms described by the cross- Wigner transform is 
well-known and considered as an asset in time-frequency analysis (e.g. radar 
theory, see Cohen [TTj , Auslander and Tolimieri [7] ) . It seems therefore that 
our approach could well open new perspectives in the topic of weak measure- 
ments and values, by importing robust techniques from these Sciences (it is 
a fact, due mainly to historical and technical reasons, that the mathematical 
techniques related to the Wigner formalism have grown faster and are more 
sophisticated in signal theory and time-frequency analysis than they are in 
quantum mechanics, so a feedback seems to be more than welcome!). 

How the weak values are related to sub-Planckian scales would also be 
interesting to investigate; the discussion in Zurek |29| . and especially the re- 
sults in Nicacio et al. [24J could certainly be useful in this context. These au- 
thors consider superpositions of an arbitrary number of Gaussian states, and 
study their motion under the action of arbitrary Hamiltonian flows. They 
show that the interference terms coming from the cross- Wigner transforms 
are always hyperbolic and survive the action of a thermal reservoir. While 
they mainly have in mind semiclassical dynamics, their approach could be 
implemented in the context of weak values. It is actually to a large extent 
sufficient to study the case of coherent states as in Subsection 12.31 because 
these states form an overcomplete set in the square-integrable functions. In 
fact, choosing an adequate lattice A of points zq in phase space the func- 
tions T(zo)£o (£,o{ x ) = (nft)~ N ^ e~\ x ^l h ~) form a Gabor frame (Grochenig 
[TTj) allowing to write an arbitrary pure state as a linear superposition of 
the states T(zq)^q. The net contribution of all cross- Wigner transforms of 
pairs (T(zq)^q,T(zi)^q) with zq ^ z\ is then the total interference leading 
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to weak values (in [29] Zurek considers a "compass state" consisting of four 
terms T(zq)S,q, of which he studies interference effects at the sub-Planckian 
scale; it would be interesting to interpret his results in terms of weak values). 

There is another aspect of the theory of weak values we have not men- 
tioned at all, if only because of lack of space and time. It is the possibility of 
reconstructing wave functions from weak values, as initiated in Lundeen et 
al. [23]. It turns out that the Wigner approach sketched in this Letter leads 
to useful formulas. For instance, on proves the following inversion formula 
(de Gosson [16], §9.4.2): Let r\ be an arbitrary square integrable function 
such that ($>\i) ^ 0; then 

2 N f 

^) = 7^TT/ WfaWTGRizoWWdzo. (28) 
W\V Jk 2N 

We can reconstruct i/j from the knowledge of the weak value provided that 
we know ((fill)- This inversion formula together with the notion of mutu- 
ally unbiased bases (MUB) could certainly play an important role in the 
reconstruction problem. 
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